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All references are textbooks and published papers 
Reference(s) used listed on each slide, lower left, as ref# 








NASA Themes and Libration Orbits 

Enterprises of Space Sciences (SSE) and Earth Sciences 
l\ are a combination of several programs and themes 
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• The Triana mission is the lone ESE mission not orbiting the Earth. 

• A major challenge is formation flying components of Constellation-X, MAXIM, TPE, and 




















Earth Science Launches 
Low Earth Orbit Formations 




- Parasol (04) 

- OCO (tbd) 
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SPECS (Submillimeter Probe of the Evolution of Cosmic Structure): 
Interferometer 1 km at L2 



Future Mission Challenges 
Considering science and operations 
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> Solar sail applications 

> Continuous control to reference trajectories 

> Tethered missions 

> Human exploration 
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FORCES ON 
PROPAGATED ORBIT 
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Gaussian Lagrange Planetary Equations 
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Potential Accelerations 
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Contains lag heating terms 
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Solar Radiation Pressure Acceleration 
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Other Perturbations 
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Inertial pointing spacecraft could have drastic changes in B c over an orbit 


Ballistic Effects 

arying the mass to area yields different decay rates 

Sample: 100kg with area of 1, 10, 25, and 50m 2 , C d =2.2 



satlOOtoI.EIapsedDays (Days) 
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Multi-Conic (patched) 

Matlab ODE 4/5 is a variable step RK 









ECEF - Navigation 
Topographic - Ground station 









8r(t) = f(r,v)~f(r , 


Describing Motion Near a Known Orbit 


<L> 



This is important since it will be our starting point for everything 
that follows 



- Hill’s - RTN 

- Clohessy-Wiltshire - RAC 

- LYLH - RIC 







Describing Motion Near a Known Orbit 
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# 5,6,7 




Transforming the EOM 












Transforming the EOM yields 
Clohessy-Wiltshire Equations 
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co . (associated with T) has twice 

cos(nt) + — sin (nt) the amplitude of the x motion 



Relative Motion 

A numerical simulation using RK8/9 and point mass 
Effect of Velocity (1 m/s) or Position(l m) Difference 



- 0.045 - 0.040 - 0.035 - 0.030 - 0.025 - 0.020 - 0.015 - 0,010 - 0.005 0.000 

main.AlongArcTrackSeparatlon (Km) 
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Graphics Ref: Collins, Meissinger, and Bell, Small Orbit Transfer Vehicle (OTV) for On-Orbit 
# 7,9 Satellite Servicing and Resupply, 15th USU Small Satellite Conference, 2001 











What Goes Wrong with an Ellipse 



Consider two trajectories r(t) and R(t). 



# 5 , 6,7 



Lambert Problem 
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‘single pass’. 






















EO-1 GSFC Formation Flying 

New Millennium Requirements 



re re 



The onboard flight software, called the EFF, shall provide the interface between the ACS / 
C&DH and the AutoCon™ system for Autonomy for transfer of all data and tables. 





















EO-1 Formation Flying Algorithm 
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where F and G are found from Gauss problem and the f _& 9 
series and C found through universal variable formulation 







EO-1 AutoCon™ Functional Description 



- Performance Objectives 
• Flight Constraints 
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Inclination Maneuver Validation : Computed AV at node crossing, of 
~ 24 cm/s (114 sec duration), Ground validation gave same results 






Difference in EO-1 Onboard 
and Ground Maneuver Three-Axis AVs 
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Formation Data from Definitive Navigation 

Solutions 
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Enables Autonomous StationKeeping, 
Formation Flying and Multiple Spacecraft Missions 
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• MAXIM Lissajous, Formation Flying of Multiple S/C 

• TPF Lissajous, Formation Flying of Multiple S/C 

(Previous missions marked in blue) 












# 1,12,13 Farquhar et al [1985] Trajectories and Orbital Maneuvers for the ISHE-3/1CE, Comet Mission, JAS 33, No. 






Goddard Space Flight Center 



350,000km, Az~ 250,000km, Class I 
Spacecraft: Mass=1254kg, Spin Stabilized, 

Notable: First Ever Multiple Gravity Assist Towards LI 










Goddard Space Flight Center 



Spacecraft: Mass=818kg, Three Axis Stabilized, 

Notable: First Gravity Assisted Constrained L2 Lissajous Orbit; Map-earth Vector 

Remains Between 0.5° and 10° off the Sun-earth Vector to Satisfy 
Communications Requirements While Avoiding Eclipses 52 
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Aspherical Configurations (Position & Rates) 
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CR3BP Analysis of Phase Space Eigenstructure 

Near Halo Orbit 







Natural Formations: 

Quasi-Periodic Relative Orbits —> 2-D Torus 









Floquet Controller 

(Remove Unstable + 2 of the 4 Center Modes) 
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Deployment into Torus 
(Remove Modes 1, 5, and 6) 
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Deployment into Natural Orbits 
(Remove Modes 1, 3, and 4) 
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Natural Formations: 
Nearly Periodic Relative Motion 
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Natural Formations: 
Expanding Vertical Orbits 



Chief S/C 








Fixed Relative Distance and Orientation 


C" 

VO 


CD 

i 

Ph 

P 

• rH 

cd 


a p 


CD 

S 

cd 

*H 

P-H 

• i-H 

t: 

CD 

a 


"Td 

CD 

X 


"id 

CD 

x 


Ph Ph 


Ph 

^H 

o 

Id 

• H 
4-> 
Ph 
CD 

S3 


n ^ 

2 td 

h—> tv j 

C 

CD 


a 

O 


xn 

C 

o 


< 

*"H +_> 

O eg 

CD Q 

CD 


CD 

O 


CD 

o 


C/5 

CD 

td 

p4 

P 

O 


C/5 

o 

Ph 

cn 

P 

o 

• H 
4—> 

cd 


cd 

$-h 

;d 

cd 

h—> 
Xfl 

• rH 

cd 

GO 

• r-H 

Ph 

P 

bD 


5o 

Q 

Q 

<P 


yp 

CD 

CD 

P 


c 

> 

> 

o 


o 

r x 

• T—1 

td 

• rH 

cd 

U 


LX 

*od 

CD 

i ■—H 

CD 

Ph 

T3 

o 

P^ 

"Td 

Id 

o 

• ^H 

5-i 

CD 

• rH 

O 

PS 

• 

f_, 

CD 

CD 

<D 

cd 

f-H 

<D 

X 

• rH 

Pg 

1 

X 

• H—1 

Ph 

O 

ph 

cd 

Ph 

Oh 

GO 

C fl 

< 

1 



Others 






















Formations Fixed in the 
Rotating Frame 



y (10 km) 
















2-S/C Formation Model 
le Sun-Earth-Moon System 



Ephemeris System = Sun+Earth+Moon Ephemeris 












Linear Targeter 








































Achievable Accuracy 
via Targeter Scheme 



Formation Distance (meters) 





Continuous Control: 

LQR vs. Input Feedback Linearization 
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Substitute x° ( t ) into EOMs and solve for u d ° ( t ) 

-+u d (t) = u;(t)-K{t){x{t)-T(t)] 

-» Apply the computed control input to the perturbed EOMs 
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Output Feedback Linearization 
(Radial Distance Control) 



Scalar Nonlinear Constraint on Control Inputs 

h(7(t),T(t)^-u (?) r r(t) = 0 





Output Feedback Linearization (OFL) 
(Radial Distance Control in the Eohemeris Model ') 
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Critically damped output response achieved in all cases 
Total AV can vary significantly for these four controllers 







OFL Control of Spherical Formations 
in the Ephemeris Model 



Relative Dynamics as Observed in the Inertial Frame 














OFL Controlled Response of Deputy S/C 
Radial Distance + Rotation Rate Tracking 
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OFL Controlled Response of Deputy S/C 
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OFL Control of Spherical Formations 
Radial Dist. + Rotation Rate 
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Inertially Fixed Formations 
in the Ephemeris Model 





Nominal Formation Keeping Cost 
(Configurations Fixed in the RLP Frame) 


















Max./Min. Cost Formations 
(Configurations Fixed in the RLP Frame) 



Nominal Relative Dynamics in the Synodic Rotating Frame 





Formation Keeping Cost Variation 
Along the SEM and L 2 Halo Families 
(Configurations Fixed in the RLP Frame) 


Q\ 

oo 




x (10 km) A (10 6 km) 

















Continuous Control in the Ephemeris Model: 
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Nearly periodic & quasi-periodic formations in the RLP frame 
Floquet controller: numerically ID solutions + stable manifolds 




Some Examples from Simulations 
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Disturbance Accommodation Model 
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Motion of Formation Flyer With Respect to 
Reference Spacecraft, in Local (S/C-l) Coordinates 
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Stellar Imager Concept 

(Using conceptual distances and control requirements 
to analyze formation possibilities) 
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Results - Libration Orbit Maintenance 
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Results - Formation Slewing 
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Results - Formation Reorientation 
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Results - Formation Reorientation 
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Results - Formation Reorientation 
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Summary 

(using example requirements and constraints) 
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comroi; ior ai couia De accomplished with better sensors to 
lessen the effect of the process and measurement noise. 


Summary 

(using example requirements and constraints) 
the controller and varying the maneuver intervals 
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• system reliability and lault detection 

• Nonlinear control and estimation 

• Second and third control tiers and new control 
strategies and algorithms. 
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DRO Mission Metrics 


co 

D 


'50 ■+- 
cd § 

6 £ 

.a = * 

r—< CT © 

— 0> y 

^ X 

^ .S T3 

O WD *3 

S «- 

v S 

1|S 

as a 


4> y 


D 

+-> 

cd 

00 

CD 

G 

X) 


ts 

(N 


o 

£ 


J- 

o 



H 

+-> 



co 



* ,-H 

T 3 

c 

o 

* t—( 

3 

£ 

1 —* 

• I —1 



o 

H 

o 

IP 

<2 . 

• 1-5 - 4 —> 

CO 

D 

J-H 

G 

§ 3 

^ CO 

V-I 

cd 

£ 

O 

CD <D 

x *- 

3 

50 

JO 

‘ 4 —> 


T3 
D 
s—» 

cd 

• l—H 

CD 

CD 

J-i 

Oh 

Oh 

cd 


co 

CD 


CD 

td 

CO 

’cd 

G 

O 

• y—t 

cd 

?-< 

(D 

Oh 

o 


Distant Retrograde Orbit (DRO) 



Asymptotic Orbit in 
the Restricted Three- 
Body Problem 



Earth Distant Retrograde Orbit (DRO) Orbit 
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DRO Formation Control Analysis 
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DRO Formation Control Analysis 
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DRO Formation Control Analysis 
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Formation Control Analysis 
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d) Resize from 25km to 50km 

e) Maintain 50km sphere (strict PD control) one month 

f) Maintain 50km sphere (loose control) one month 
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DRO Formation Control Analysis 

> Earth/Moon L4 Libration Orbit 

> Spacecraft controlled to maintain only relative separations 
^ Plots show formation position and drift (sphere represent 2 
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General Theory of Decentralized Control 



•Data transmission requirements are * W 

MINIMIZED. Node 4 
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DST/Numerical Comparisons 
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Libration Point Trajectory Generation Process 





















































irajectory ueneratec 
with Manifold States 







JWST DST Perspective 





1 raj ectory Lrenerated starting \ 

with Manifold States Starting Point 
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Coordinate system transformation 
Orbital to inertial coordinates 
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Principles Behind Decentralized Control 










The LQG Decentralized Controller Overview 



TO NETWORK 





